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ABSTRACT 

There is some evidence, though yet unconfirmed, that Sagittarius A* — the 
supermassive black hole at the Galactic center — emits its radio waves modu- 
lated with a ~ 100-day period. What is intriguing about this apparent quasi- 
periodicity is that, though the amplitude of the modulation increases with de- 
creasing wavelength (from 3.6 to 1.3 cm), the quasi-period itself does not seem to 
depend on the frequency of the radiation. It is difficult to imagine how a binary 
companion, were that the cause of this modulation, could have escaped detec- 
tion until now. Instead, it has been suggested that the spin-induced precession 
of a disk surrounding a slowly rotating black hole could have the right period 
to account for this behavior. In this paper, we examine how Sagittarius A*'s 
light curve could be modulated by this mechanism. We demonstrate that the 
partial occultation of a nonthermal halo by a compact, radio-opaque disk does 
indeed produce the observed frequency- dependent amplitude. This appears to 
be in line with other observational arguments suggesting that Sagittarius A*'s 
mm/sub-mm spectrum is produced by a ~ 10 Schwarzschild-radius disk, whereas 
its cm-waves originate from a nonthermal particle distribution in a halo extend- 
ing out to over 20 Schwarzschild radii. Interestingly, this model suggests that 
the observed period corresponds to half the precession period and that a non- 
axisymmetric disk could produce a second period roughly twice as long as the 
first. 



Subject headings: accretion, accretion disks - - black hole physics - - Galaxy: 
center — gravitation — radiation mechanisms: nonthermal — relativity 
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1. Introduction 

The weekly VLA monitoring of Sagittarius A*, the ~ 4x 10 6 M supermassive black hole 
at the center of our Galaxy, has accumulated over 20 years of variability data at 1.3, 2.0, 3.6, 
6.0, and 20 cm wavelengths. The sampling within this period has been somewhat irregular 
(Zhao, Bower, and Goss 2001). Nonetheless, the power spectral density (PSD) reveals a 
clear peak near 1 x 10 -7 Hz, with a progressively smaller significance at longer wavelengths. 
This frequency corresponds to a periodic modulation of approximately 100 — 120 days; the 
actual best-fit period extracted from the combined data sets is 106 days. 

This result is at once intriguing and unsettling. The fact that a ~ 17-minute Keplerian 
period has been seen in Sagittarius A*'s infrared emission (Genzel et al. 2003) makes this 
cyclic modulation easier to accept. Yet the implied radio period (~ 106 days) contrasts 
sharply with the dynamical time scale (~ 17 — 20 minutes) associated with motion in the 
inner disk (see Melia and Falcke 2001 for a recent review). 

Perhaps the periodic radio signal is a false detection due to a combination of a random 
process and the irregular sampling pattern. However Monte Carlo tests with data created 
from various sources of noise using this same sampling don't seem to bear this out. Regardless 
of the type of noise used in the simulations — including white noise, Gaussian noise around a 
mean, and a Poisson distribution of flares — the probability of false detection due to any such 
random process appears to be less than 5%. (However, one should keep in mind that 1 / / — or 
red — noise could in principle mimic such a periodic signal if the random fluctuations have 
an appropriate scale; G. Bower, private communication.) Perhaps also supporting the view 
that this period may be real is the observation that both the absolute (AS') and fractional 
(AS/S) amplitudes of the pulsed component increase toward shorter wavelengths, yet the 
period appears to be independent of wavelength. 

The 106-day cycle evident primarily at 1.3 and 2.0 cm may be a valuable tool for probing 
Sagittarius A*'s inner workings should it truly have something to do with the source. High- 
resolution VLA observations have already ruled out the possibility that such a period might 
be produced by an orbiting emitting object (Bower and Backer 1998). The 106-day orbit 
of a companion to Sagittarius A* would have a radius ~ 60 A.U., corresponding to an 
angular separation of ~ 8 milliarcseconds at 8 kpc. A compact 0.2-Jy source separated 
from Sagittarius A* by this amount would have easily been detected with the VLBA at 
wavelengths shorter than 3.6 cm. The unlikelihood of Sagittarius A* having an orbiting 
companion is further supported by its observed lack of proper motion, which precludes any 
possible association with rapidly moving components. In addition, a stellar origin for such 
a source would fall well short of the power required to account for the measured radio 
emission. All in all, the evidence seems to favor an interpretation in which the 106-day 
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periodic variations, if real, are intrinsic to Sagittarius A* itself. 

The characteristics of this 106-day cycle constrain the nature of its origin rather tightly. 
First, the observed period is, as we have said, independent of wavelength. The emission 
in Sagittarius A* at different frequencies is produced on different spatial scales (see, e.g., 
Melia, Jokipii, and Narayanan 1992), so the period should be induced by a single process. 
Otherwise, we would expect to see different periods at different frequencies. What is required 
is something that can cause correlated fluctuations across a broad range of wavelengths. 
Second, the period is four orders of magnitude longer than the dynamical time scale in 
the inner disk surrounding Sagittarius A*. Could it be produced on a much larger spatial 
scale? Sagittarius A*'s 2-cm emission is produced within a region no bigger than ~ 24 A.U. 
(Krichbaum et al. 1999), for which the corresponding dynamical time scale is calculated to 
be about one and a half days; so the answer is apparently no. Higher frequency emission is 
produced within still smaller regions, associated with even smaller time scales. 

We may ask then, whether this modulation could be produced by a corrugation wave 
in an accretion disk, which is used to account for the quasi-periodic oscillations (QPOs) 
seen in low-mass X-ray binaries. These waves have periods that are much longer than the 
corresponding dynamical time scale (Kato 1990), but they depend on radius and thus may 
not be able to account for the first feature described above. Moreover, Sagittarius A*'s 
light curves show quite stable periodic fluctuations, rather than the uncorrelated segments 
constituting QPOs. The evidence is pointing to a single process evolving in a relatively 
confined region, certainly no bigger than 100 Schwarzschild radii (r$ = IGMjc 2 « 10 12 cm 
for a 4 x 10 6 M black hole). 

In an earlier paper (Liu and Melia 2002), we noted that the gravitational acceleration in a 
Kerr metric acquires a dependence on poloidal angle (relative to the black hole's spin axis), so 
that matter orbiting above or below the equatorial plane experiences a restoring force toward 
the equator. This results in the precession of its angular momentum vector about the black 
hole's axis of rotation. The physical conditions in Sagittarius A*'s Keplerian region lead to 
strong coupling between neighboring rings, and the ~ 20 — 30 r s compact disk therefore 
torques more or less as a rigid body. Under an appropriate set of circumstances (see Liu 
and Melia 2002), the precession period can exceed a hundred days, and the long-term radio 
modulation in Sagittarius A* may be closely related to its short-term X-ray and infrared 
variability after all — via the dynamical properties of the disk in a rotating spacetime. 

However, whereas this earlier treatment established the viability of a spin-induced disk 
precession in accounting for the 106-day modulation, it left open the question of how exactly 
this periodicity is manifested. The purpose of this paper is to demonstrate that the observed 
fluctuations, frequency- dependent amplitudes, and periodicity can in fact be produced as a 
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result of partial occultation of a nonthermal halo surrounding Sagittarius A* by the pivoting 
disk. As we shall see, the properties inferred for this source based on its spectrum and 
polarization characteristics (see Liu and Melia 2001), produce a surprisingly accurate fit to 
the observed radio lightcurve, its modulated amplitude, and the frequency-dependent signal 
strength. In § 2 we describe the method used for this analysis, including the source geometry 
and particle properties. We summarize our results in § 3, and present our conclusions in § 4. 

2. The Method 

The level of polarization seen in Sagittarius A* at mm/sub-mm wavelengths approaches 
7—10% (Aitken et al. 2000). However, this object reveals a lack (< 1%) of linear polarization 
below 112 GHz, though some circular polarization (~ 1%) has been detected (Bower et al. 
1999; Bower et al. 2001). These prominent spectral and polarimetric differences (Melia, 
Bromley, & Liu 2001) between the cm and the mm/sub-mm bands suggest two different 
emission components in Sgr A*. As we have already noted, higher frequencies correspond 
to smaller spatial scales (see also Narayan et al. 1995), so the mm/sub-mm radiation is 
likely produced in the vicinity of the black hole. Earlier work (e.g., Melia 1992, 1994; Coker 
& Melia 1997) has indicated that Sgr A* is accreting from the stellar winds surrounding 
the black hole and that the infalling gas circularizes at a radius of ~ 20 — 800 r$- Recent 
work on Sagittarius A*'s emissivity (Melia, Liu, & Coker 2001; Bromley et al. 2001) has 
demonstrated that the inner 10 r$ of the resultant Keplerian structure can not only account 
for the mm/sub-mm properties via thermal synchrotron emission, but it may also produce 
Sgr A*'s X-ray spectrum in the quiescent state (Baganoff et al. 2001) via Comptonization of 
the mm/sub-mm photons. On the other hand, the cm radio emission appears to be produced 
by non-thermal synchrotron emission (Liu & Melia 2001). 

The integrated ~ 1 — 20 cm luminosity of Sgr A* is comparable to the power extracted 
from its spin energy via a Blandford-Znajek type of electromagnetic process if a/M 0.1 
(Liu and Melia 2002). The 106-day modulation could therefore presumably arise when the 
precessing disk periodically shadows the non-thermal particles flooding the region surround- 
ing the black hole as they escape from their creation site near the event horizon, essentially 
forming an expanding halo of relativistic particles. Some observational evidence for this has 
recently been provided by VLBA closure amplitude imaging techniques at 7 mm (Bower 
et al. 2004), which point to an intrinsic source size of ~ 25 r s for Sagittarius A* at this 
wavelength. 

We will therefore adopt the basic model displayed in Figure 1. The compact disk (with 
outer radius Rd) that produces the mm/sub-mm emission is opaque to radiation longward 
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of ~ 1 cm, which we assume originates from the surrounding, semi-transparent halo with 
radius Rh- To accurately determine the disk shadowing effect, we integrate the non-thermal 
synchrotron emissivity j v (erg cm -3 s -1 ster -1 Hz -1 ) along each given line-of-sight (see Figure 
1), including the effects of opacity, such that 

Iu{s)= [' M<r) e-r-W dff . (1) 

J so 

In this expression, 

t u (s) = / a v {a) da (2) 

J SO 

is the optical depth in terms of the absorption coefficient a u (a) (cm -1 ). The flux is then 
calculated by integrating l v over all solid angles pertaining to the source: 



? = J I v cos 9 dQ, (3) 



where 9 is the angle between the line of sight to the black hole (at the center of this geometric 
configuration) and the emitting area associated with dQ. 



Numerically, we calculate the flux on a grid (see Figure 2), assuming F u to be constant 
over each mesh element, so that the total flux may be written 

d 2 



where d is the distance to the Galactic center, e is the side length of each mesh element, and 
9 « for e/d < I. 



2.1. Power Law Distribution 

We assume that within the halo the emitting particles form a power-law distribution 

N(E) dE = N E~ p dE , (5) 

with p — 2. Thus, 

\/3 e 3 r°° 

jv — 7 ^.Bsin^ / N{E) F{x) dE , (6) 

47r m& J 

where B is the magnetic field strength, ip is the pitch angle (the angle between the electron's 
velocity and B), F(x) = x if 5 / 3 (£) d^, and K 5 /z is the modified Bessel function of the 
second kind. Similarly, 



87rm v 2 r J dE \ E 2 J 
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Equations (6) and (7) then become 



a 



3v = 



V 



dN {B sin^)**' - * ? 
C 2 N {B sin^) 2 //" 3 , 



(8) 
(9) 



where 



Ci = 




2m 5 c 9 ' 



(10) 



C 2 = 



67r m 4 c 5 



(11) 



Equations (8) and (9) are used in Equations (1) through (4) to calculate the total flux. 



For the calculations described here, the disk precession is handled analytically, with an 
orientation prescribed as a function of time, based on our previous work (Liu and Melia 2002). 
A full SPH simulation is currently underway, and a more realistic time-dependent geometric 
profile will be published elsewhere (Rockefeller, Fryer, and Melia 2005). A key concern is 
whether the disk succumbs to the Bardeen-Pettersen effect, in which the radially-dependent 
precession frequency can lead to dissipation between neighboring rings, ultimately forcing 
the inner ~ 10 — 20 r$ of a disk inclined to the black hole's spin axis to eventually settle into 
the equatorial plane. It is commonly thought that an accretion disk surrounding a spinning 
black hole must be warped, with an overall angular momentum vector possibly misaligned 
relative to the spin axis, but with its inner portion fully flattened at the equator. But detailed 
hydrodynamical simulations are now showing that this effect, although fully realized under 
a majority of physical conditions, can be negated in cases where other couplings between 
neighboring rings in the disk are strong. An example of this occurs when the disk plasma has 
a small Mach number (typically < 5), for which large pressure gradients can then distribute 
the radially- dependent precession torque and force the disk to rotate as a rigid body (see, 
e.g., Nelson & Papaloizou 2000). In the case of Sgr A*, we estimate from the spectrum 
and polarization properties of its disk that the Mach number in this system is ~ 3. Not 
surprisingly then, the latest SPH simulations confirm the earlier semi-analytic results that 
the small disk in Sgr A* is probably precessing as a rigid body. In our calculation, we 
therefore evolve the disk orientation analytically in time, assuming rigid body precession, 
and recalculate the flux at 1.3, 2.0, and 3.6 cm for each step. 



2.2. Parameters and model assumptions 



As we can see from Equations (8) and (9), the model contains five basic parameters: N , 
B, R h , R d and 7. Here 7 is the angle between the angular momentum vector of the disc and 
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the spin axis of the black hole (see Fig. 1). The other two angles a and (3 change in every 
time step. The pitch angle ip depends on distance a in Equation (1). For relativistic electron 
energies, we can approximate ip ~ i?, where $ is the angle between the magnetic field vector 
and the line of sight. For the sake of simplicity, iVo is taken to be a power-law function oc r~ £ 
(with e — 0, 1, 2), in terms of r, the distance from the center of the halo. The magnetic field 
B is assumed constant throughout the halo. Consequently, there are effectively six overall 
parameters in this model: N (r = Rh), B, and e, Rh, Rj and 7. 

Throughout this work the accretion disk is considered to be optically thick to all radia- 
tion onward of 1 cm, so any ray of light intercepting the disk is immediately stopped at that 
point. Since the distance to the Galactic center d is very large compared to the radius of the 
halo we may also assume parallel lines of sight. 

We have explored models with R^ ~ 17 r$, Rh ~ 22 for the halo, and a tilting angle 
7 between L and S of ~ 60° (or 30° measured from the xy-plane). Based on earlier spectral 
fitting calculations, we know that typical values for N (r = R h ) are of order 10 14 cm -3 erg p . 
For a given power-law dependence of N , this then also fixes the density throughout the 
emitting region. Also based on earlier spectral fitting calculations (see Liu and Melia 2001), 
we infer a magnetic field intensity of order 1 G. 

We point out, however, that the specific choice of the function N (r) has little effect on 
the shape of the light curve (as we shall see in Figs. 3 to 8 below), though the calculated 
flux does change in response to the changing column depth through the emitting medium. 
Consequently, we have chosen the model assumptions to be as simple as possible. A uniform 
electron distribution with e = 0, though unrealistic, allows a first critical analysis of the 
shadowing effect. 

3. Results 

A good fit to the folded light curves (shown in Figs. 3 and 4) is possible with either a 
uniform electron distribution, or a power law. For the sake of specificity, we here focus on the 
results produced with a uniform halo (Figs. 3, 4, and 5), and one in which Nq(t) oc 1/r (Figs. 
6, 7, and 8). We have also found that the light curve produced as a result of shadowing (dark 
solid curve in these figures) is only weakly dependent on the absolute value of N Q (Rh) and 
B, which instead affect primarily the absolute flux level. The shadowing effect is therefore 
relatively independent of the details of the underlying model; it appears to be a rather robust 
phenomenon. 

In Figure 3, we show the calculated light curve in comparison with the data at 1.3 cm. 
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This fit was produced with N (R h ) = 3.6 x 10 13 cm" 3 erg 2 , B = 0.8 G, R h = 21.5 r s and 
Rd = 17.5 rs- Every data point on the solid line represents the flux (in Jy) obtained for a 
certain position of the accretion disk, which precesses around the (fixed) spin vector of the 
black hole. In every time step between time days (which also corresponds to 0°) and 106 
days (which corresponds to an angle of 180° of L with respect to the starting point) the 
precessing angle is increased. At the maximum 84 days), the accretion disk is seen edge- 
on and consequently, since the disk is here assumed to have no thickness, the integration 
over the synchrotron emissivity extends over the full extent of the halo. In order to illustrate 
this, for a tilting angle of 7 = 90 degrees, i.e., L parallel to the x-axis in Figure 1, the curve 
produced by our simulation would be constant at the value of the maximum in Figure 3, 
since the accretion disk would be seen edge-on in every time step. 

It should be emphasized that in Figure 3, the mean flux has not been adjusted vertically, 
but rather is set by the parameter values suggested by earlier spectral studies. One of the 
principal results of our simulation is that the parameters that yield this correct flux, together 
with a reasonable tilting angle of 60°, yield both the correct amplitude of the observed 
modulation, and its variation in time. 

In Figure 3, we also show the corresponding light curve for 2.0 cm, using the same 
parameters as those for 1.3 cm. Notice that here the shadowing is less effective, which is due 
to two factors: first, the halo becomes progressively more opaque with increasing wavelength, 
which means that a greater fraction of the observed intensity arises from the medium in front 
of the disk; second, the emitting volume contributing to the overall flux increases, rendering 
the fraction of the halo occulted by the disk less significant with increasing wavelength. The 
optical depth through the halo at 1.3 cm is r v 10 _1 — 10~ 2 , so most of the photons pass 
through the medium unaffected. However, r v goes as z/~ 3 (see Eq. 9); at 2.0 cm, it is > 10 _1 . 
We note that both the shape of the light curve and the amplitude of modulation fit the data 
just as well at 2.0 cm as they do at 1.3 cm, for the same system parameters. 

By 3.6 cm, r„ has increased to values > 1 and, not surprisingly, most of the modulation 
due to the precessing disk, which lies below the photosphere at this wavelength, has disap- 
peared (see Fig. 5). This trend is reflected in the data as well, though we note that our 
calculated flux is a factor 2 lower than that observed. We attribute this to the fact that for 
simplicity we have assumed a fixed halo size. In reality, the fact that t v is increasing with 
increasing wavelength also means that the emitting volume contributing to the overall flux 
must itself increase. It would be a simple matter to reproduce the correct flux at 3.6 cm 
by making R h wavelength dependent, but this would violate our goal of keeping this model 
as simple as possible. In principle, a full numerical simulation of the disk precession and 
particle acceleration and escape should yield the correct behavior of quantities such as iVo 
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with radius, which would then make this point moot. 

Figure 4 demonstrates the same effects as those discussed in Figure 3, except now for 
a halo with uniform magnetic field, but with a density N (r) oc 1/r, with N (R h ) = 2.45 x 
10 13 cm~ 3 erg 2 and with B = 0.8 G. The fact that a reasonable fit to the data may be made 
with such disparate halo geometries affirms the robustness of this model. However, changing 
the parameter values does produce quantitative variations in the modulation amplitude, and 
its dependence on wavelength. Figure 5 shows the sensitivity of our simulation to changes 
in No(Rh). We see that the shape of the curve (at 1.3 cm) is not affected. This is hardly 
surprising, since r v < 1 for all these densities. However, since j u oc N , the overall flux does 
change with increasing iVo- 

The consequence of changing the magnetic field is shown in Figure 5. As was the case 
for N , the shape of the curve does not depend on this value. Also, since the dependence 
of a u and j u on B is only to the power 3/2, the impact on the average flux is also less 
pronounced. 

The impact of parameter variations on the light curve at other frequencies is similar to 
that for 1.3 cm, so we won't show these here. 

Changes in the disk's outer radius produce the modified modulations shown in Figure 
5. All the curves coincide at the peak since the edge-on disk does not shadow any part of 
the halo. Also, the amplitude decreases as the disk becomes smaller. 

4. Conclusions 

The picture described here works rather well in accounting for (1) the absolute flux, (2) 
the observed modulation, and (3) the amplitude of the observed peaks in the 1.3 and 2.0 cm 
light curves of Sagittarius A*. However, future work should include a more realistic time- 
dependent disk profile, the halo, and a more detailed exploration of possible magnetic field 
configurations. For example, in this paper, we have restricted our attention exclusively to 
isotropic halos surrounding Sgr A*, though we have allowed for a possible radial variation in 
its physical conditions, such as the particle number density and the magnetic field. Clearly, 
though, one might expect some impact on the results should the halo not be spherical, but 
rather toroidal, or cylindrical. But as we have already indicated earlier, the halo is optically 
thin to radio emission, and the dominant effect is the shadowing of this emitter by the 
disk. Our brief survey of the various possible geometries has indicated that the shape of the 
modulated light curve is only weakly dependent on the halo's structure. Nonetheless, the 
absolute flux level is indeed sensitive to the halo's geometry. Some of the additional work 
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necessary to fully explore the range of possibilities (e.g., with SPH simulations) is currently 
underway and will be reported elsewhere. 

Finally, we point out here something that was overlooked in Liu and Melia (2002). This 
is the fact that if the disk precession is responsible for the modulation, then the observed 
period should actually correspond to half of the precession period because of symmetry. This 
has some impact on the inferred value of the black hole spin 2a/ Vg. In addition, it may be 
possible to see two periods, the primary one associated with half the precession period, and 
a second at roughly twice the first, corresponding to the full period. The latter would arise 
if the disk were not perfectly symmetric, so that the tilt in the first half of a cycle is not 
exactly the same as the tilt in the second. This is a prediction of the model presented in 
this paper and could be included in future work. Interestingly, there may already be some 
evidence for a second period, roughly twice as long as the first, in data acquired since 2000 
(Zhao, private communication). 

To close, it is still not entirely clear whether this ~ 100-day period in Sagittarius A*'s 
emission is in fact due to an intrinsic modulation. However, we have now shown that a 
spin-induced disk precession for a slowly rotating black hole can account for this long period 
and, in addition, that the partial occultation of a ~ 1 — 3.6-cm emitting halo by the radio- 
opaque precessing disk can also produce the correct amplitude and time-dependent light 
curves seen at 1.3 and 2.0 cm. The halo becomes optically thick longward of ~ 2.0 cm, so the 
shadowing by the disk can no longer produce a modulation at these wavelengths. Future work 
observationally should focus more carefully on the question of whether red noise fluctuations 
could mimic the effect of a long periodic modulation in this source. Theoretically, future 
work will determine the temporal profile of the precessing disk more accurately, and will 
provide us with a better grasp of the halo's internal structure. 

M.P. would like to thank the German Fulbright Commission for support far beyond the 
financial aspect of academics, which made this work possible. Furthermore he likes to thank 
Steward Observatory for kind hospitality and Casey Meakin and Philipp Strack for helpful 
discussions. This work was supported by NASA grant NAG5-9205 and NSF grant 0402502 
at the University of Arizona. 
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Fig. 1. — The geometry used in this model. Here, d denotes the distance to the Galactic 
center, R h is the radius of the nonthermal halo surrounding Sgr A*, R d the radius of the 
compact accretion disk, S is the spin vector of the black hole, and L is the angular momentum 
vector of the disk making angles a, (3, and 7 to the x, y, z axes, respectively. 
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Fig. 3. — Top panel: calculated light curve (solid) at 1.3 cm, assuming a uniform nonthermal 
halo with radius Rh = 21.5 r s , enclosing a compact disk with radius Rd = 17.5 r S - See text 
for additional parameter values. The data are folded over a 106-day period, and are taken 
from Zhao et al. (2001). Middle panel: same as above, except here at a wavelength of 2.0 
cm. Bottom panel: same as above, except here at a wavelength of 3.6 cm. 
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Fig. 4. — Same as Figure 3, except here for a halo with a 1/r density profile. See text for 
additional parameter values. 
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Fig. 5. — Top panel: the effect on Sagittarius A*'s light curve at f.3 cm due to variations 
(as shown) in the halo's nonthermal particle density. Middle panel: the effect on Sagittarius 
A*'s light curve at f.3 cm due to variations (as shown) in the halo's magnetic field. Bottom 
panel: the effect on Sagittarius A*'s light curve at f.3 cm due to variations (as shown) in 
the disk size. 



